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Abstract 

The relative entropy in two-dimensional field theory is studied on a cylin- 
der geometry, interpreted as finite-temperature field theory. The width of the 
■ cylinder provides an infrared scale that allows us to define a dimensionless rel- 

ative entropy analogous to Zamolodchikov's c function. The one-dimensional 
\ quantum thermodynamic entropy gives rise to another monotonic dimensionless 

i-^ ■ quantity. I illustrate these monotonicity theorems with examples ranging from 

£> . free field theories to interacting models soluble with the thermodynamic Bethe 

ansatz. Both dimensionless entropies are explicitly shown to be monotonic in 
$h ' the examples that we analyze. 
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1 Introduction 

It has been shown that the irreversible character of the renormalization group (RG) 
can be cast in a sort of H theorem analogous to Boltzmann's, thus generalizing this 
theorem from ordinary time evolution to the evolution with the RG parameter Q. 
The irreversible quantity, the field theory entropy relative to a fixed point of the RG, 
is a monotonic function of the coupling constants and increases in the crossover from 
one fixed point to another less stable. However, the Wilson RG picture considered 
in U, wherein one has to deal with all the couplings generated by the RG action, 
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relevant and irrelevant alike, turns out to be too complex and was indeed assimilated 
to a non-equilibrium thermodynamics setting. One can start with only the relevant 
couplings but then one must utilize a different RG which changes some infrared (IR) 
scale. A possibility is to define the field theory on a finite geometry characterized by 
some parameter, loosely associated to its size, which plays the role of IR scale. Then 
the monotonicity theorem for the relative entropy can be cast as a RG theorem similar 
to the celebrated Zamolodchikov c theorem || . 

Among the various geometries we could consider, the cylinder stands out for its 
simplicity. It is defined by only one scale, the length of the compact dimension, and 
the finite-size corrections to the partition function turn out to be computable. More- 
over, on a cylinder of circumference j3 the monotonicity theorem adopts a form with a 
thermodynamic interpretation, the temperature being T = 1/(3 M. Thus the inverse 
temperature is used as RG parameter, providing a thermodynamic interpretation of 
the RG, as in Ref. |§. The connection with concepts of 1 + 1 QFT at finite temper- 
ature is intellectually appealing and useful for computational purposes. For example, 
finite-size corrections are calculated in terms of the properties of one- dimensional (Id) 
quantum gases. In addition to the relative entropy, the Id quantum entropy provides 
another monotonic quantity with a different interpretation. We must remark that the 
definition and monotonicity of the relative entropy, as exposed in Ref. [|IJ, already 
have a thermodynamical motivation in the 2d context, independently of the type of 
geometry. However, in field theory we prefer to dissociate the coupling from a thermal 
interpretation and we reserve the concept of temperature for its role in the Id quantum 
picture. Nevertheless, it shall be evident that the proofs of the monotonicity theorems 
for the 2d relative entropy or for the Id quantum entropy are essentially the same. 

In Ref. [§] these ideas were illustrated only with free field models and the calcu- 
lations of the corresponding finite-size corrections were presented very concisely. We 
shall begin here with a more detailed analysis of the properties of both types of entropy, 
in particular, considering whether they are universal quantities. Next, we proceed to 
the explicit calculation of the finite-size corrections for soluble models corresponding to 
free-field theories, including thermodynamic quantities as well as the expectation values 
of the stress tensor, and hence of the entropic monotonic quantities. The properties of 
these quantities will be displayed in the corresponding plottings. Further to free-field 
models, it will be demonstrated that interacting models are also suitable for calcula- 
tion of their finite-size corrections and monotonic quantities with powerful methods. 
In particular, integrable models on the cylinder are appropriate for application of the 
thermodynamic Bethe ansatz (TBA). Plots of the monotonic quantities obtained with 
this method display similar behaviour to those of free-field models. 

The paper is divided in three parts. The first part is devoted to formulating the 
monotonicity theorems for 2c? field theory and to giving its thermodynamic interpreta- 
tion on the cylinder. The second part applies these theorems to the relatively simple 
cases of the Gaussian and Ising models. They allow for an explicit calculation of ther- 
modynamic quantities and their connection with the components of the stress tensor. 
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Section 3 is devoted to interacting models which lend themselves to computation of 
thermodynamic quantities. The essential tool is the thermodynamic Bethe ansatz, 
which is first applied to models with purely statistical interaction, resulting again in 
explicit expressions for the relevant quantities, and in second place to models in which 
the TBA equations have to be solved numerically. Afterwards there comes a discussion 
of the results obtained and, finally, two appendices, the first one on the method for 
the computation of finite-size corrections based on the Euler-MacLaurin formula and 
the second one on the computation of the expectation value of the stress tensor on the 
cylinder for free theories. 



2 Entropic C theorems 

2.1 General properties of the relative entropy in two-dimensional 
field theory 

Let us briefly recall some concepts already introduced in Ref. [p]]. The field theory 
probability distribution associated to some statistical system is given by 



V[(f),{X}] = e 



-7[</>,{A}]+W[{A}] 



and depends on some stochastic field <fi and a set of coupling constants {A}. The 
quantity ^[{A}] is needed for normalization and is of course minus the logarithm of 
the partition function. A composite field is defined as the derivative of the action with 
respect to some coupling constant, 

h - |. (2) 

For example, if we consider the thermal coupling, the coupling constant is the inverse 
temperature and the composite field represents the energy. As is usual, we assume for 
simplicity that I [(f), {A}] is linear in the coupling constants. 

The relative entropy, a concept borrowed from probability theory, turns out to be 
the Legendre transform of ^(A) — 1^(0) with respect to A Q: 

S Iel {X) = W(\) - W(0) - X^- = W-W -\ (fx) ■ (3) 

dX 

Obviously, S^O) = 0. Furthermore, as a straightforward consequence of its definition, 

dS rel dW d fdW 

A— ^ = A— -A— A 



dX dX dX \ dX 
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For the thermal coupling, S re \ has indeed the interpretation of a real thermodynamic 
entropy which increases with temperature. In other cases, it may or may not have a 
thermodynamic interpretation but its properties hold nonetheless. 

Some qualifications are in order. In field theory we deal with local fields, so f\ = 
J $a, where &\(z) is a local composite field, function of the 2d coordinates z = x% + 
ix 2 - We must remark that, although these fields are usually constructed as actual 
composites of the basic field 0, the existence of this field needs not be assumed, as 
in some modern formulations where it is replaced by the action principle ||. This 
remark is important when we start from a 2d conformal field theory. To prevent the 
appearance of ultraviolet (UV) divergences it is convenient to define all the quantities 
with a UV cutoff A (for example, W[X, A]) which must be eventually removed to define 
universal quantities. Even though W is non universal, we expect S Te \ to be . In order 
to have universality, we consider RG relevant or marginal couplings: In two dimensions 
the scaling dimension of the field $ must be such that < dq> < 2. This condition 
may not be sufficient and shall be made more precise shortly. W and S re \ are extensive 
and it is convenient to define the associated specific quantities dividing by the total 
volume — or area in two dimensions. Henceforth, we use specific quantities but keep 
the same notation for simplicity. We are interested in an entropy relative to a RG fixed 
point, so we must substract from the coupling constants their values at that point. 
(The fixed-point coupling constants may be null in some cases.) Finally, there is an 
assumption of positivity of the probability distribution implied in the inequality (|]), 
like in Zamolodchikov's theorem. 

To derive a universal expression for the specific S Te \ we must analyse its dependence 
on the UV cutoff. We can use the scaling form of the specific W, 

W/(A,A) = A 2 ^— J, (5) 

where y = 2 — d® > is the dimension of the coupling A. For the thermal field, 
the local energy density, y is the inverse of the critical exponent v. If the scaling 
function is continuously differentiable around zero (class C 1 ), and we denote F Q = J^(0), 
Fi = ^'(0), W can be expanded as 

W(X, A) = A 2 F + Fx X 2/y + A 2 o(A~ 2 ), (6) 

with o(A~ 2 ) asymptotically smaller than A -2 , hence resulting in a vanishing term as 
A — > oo. Given that the UV divergent term of this expansion cancels in W^(A) — 1^(0), 
the relative entropy yields a finite result in the infinite cutoff limit, namely, 

sua) = w(x) - w(o) -x d ^ = F l y — x 2 'y. (7) 

dX y 

Thus the significance of the assumed regularity condition on the scaling function is 
that it is sufficient to endow the monotonicity theorem with universality. One can 
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certainly think of simple functions that are not class C . For example, the function 
T{x) = F Q — x lnx + o(x), which will appear in some of the models studied later. 

We now examine the question of universality in terms of local fields. This method 
will lead us to a more concrete formulation. Let us begin by writing the monotonicity 
theorem (gj) as 

A^ = A 2 |d 2 z(:$(z)::$(0):)>0, (8) 

with the use of the definition of normal-ordered composite fields, :$: = $— ($). We 
can study the UV convergence of this integral. As a prerequisite, note that possible UV 
divergences in the definition of the composite field $ are removed by the substraction 
of ($). The most singular part of the correlation function for short distance is given 
by 

(:$(z)::$(0):) ~ |z|~ 2d *. (9) 

Hence, the integral converges if < d$ < 1, that is, 1 < y < 2. Then the derivative 
of the relative entropy, dS re i/d\, is a universal quantity and so is S re \, because the 
integration constant is fixed by the condition S^O) = 0. For dimensional reasons, it 
must adopt a form like that in (^): 

s vcl (\) = B\ 2 /y, (io) 

where B is a constant. In fact, upon inversion of the Legendre transform this form 
implies that W has the previous first-order expansion except in the case of y — 2. 

Fields $ satisfying < < 1 are called strongly relevant || |6|. They include 
the thermal coupling of the unitary minimal models of conformal field theory (CFT), 
except the Ising model, wherein the local energy density has d$ = 1. We shall find that 
the relative entropy of the Ising model is indeed non universal. In principle, fields with 
1 < d$ < 2 give rise to a non-universal relative entropy, S re \(X, A). It is monotonic and 
essentially independent of A as long as X 2 ^ y <C A, which is the condition necessary for 
the continuum field theory of the statistical system to be meaningful. In this sense, 
one may consider this non-universal relative entropy within the philosophy of effective 
field theories, a term which refers to theories that are not renormalizable but suitable 
for calculation of many physical quantities for scales much lower than the cutoff. 

We must remark that the simple power-law forms of the relative entropy (|T0D and 
the monotonicity theorem are not very informative, in the sense that, once we know 
that S re \ is finite, they follow from dimensional analysis. We thus see the necessity 
of introducing a new parameter, for example, through a finite geometry. We will 
indeed obtain a richer and more illuminating version of the relative entropy and the 
monotonicity theorem when we introduce a finite geometry. 

Let us introduce the stress tensor trace, 9 := T£. Since 6 gives the response 
to a change of scale and the only scale is in the coupling constant, it is in general 
proportional to the relevant field $: 

9 = yA$. 



5 



Hence, we can put the monotonicity theorem for the specific relative entropy in an 
interesting form: 

A ^L(A,A) = 1 J#z{:e(z)::Q{0):)>0. (11) 
We will have the occasion to comment on this form in what follows. 

2.2 Finite-size corrections. The cylinder and one-dimensional 
thermodynamics. 

So far, the relative entropy has been proved to be monotonic with respect to the 
coupling constants. Now we would like to reformulate the monotonicity theorem for 
the relative entropy as showing irreversibility under the RG. We need to substitute the 
coupling constant A by some quantity which can be interpreted as a RG parameter. 
A common way to introduce a RG parameter is through some IR scale. For example, 
we may consider a finite size system with a characteristic length, such as a strip or 
cylinder of width L. According to finite-size scaling ideas, the free energy can be split 
into a bulk part and a universal finite-size correction. Tha latter constitutes a suitable 
function to derive a non-trivial relative entropy. Moreover, one can take advantage of 
the fact that the classical partition function on a cylinder of width (3 is equivalent to 
the one-dimensional quantum partition function at temperature T — 1/(3 to give the 
RG a thermodynamic interpretation [RJ. Indeed, relevant thermodynamic functions of 
this quantum system are given by derivatives with respect to (3. The first one is the 
energy, which has one part independent of the temperature and another that vanishes 
at T = 0, corresponding to the bulk part and the finite size correction, respectively. The 
part independent of the temperature, which is non universal, represents the ground- 
state energy. A more interesting quantity is the specific one-dimensional quantum 
entropy, which turns out to be universal and will prove to be the right quantity for a 
thermodynamic monotonicity theorem. 

Let us then consider the system on a cylinder, equivalent to finite temperature 
field theory. The partition function is Z = Tre~@ H , which can be represented as a 
functional integral on S 1 x IR with (3 = 1/T the length of the compact dimension. We 
assume that the specific logarithm of the partition function on a cylinder of width (3 
and length L as L — > oo can be split into a bulk part and a finite-size correction^ 

zM = 4 = eo(A , A) , + £(M, (12) 

where C((3,X) is a universal dimensionless function having a finite limit as (3 — > oo. 
Hence, defining x = (3 X 1 ^ we write C(/3,A) as a single- variable function, C(x). At a 
RG fixed point it is proportional to the CFT central charge, C(0) = — nc/6 J7|, §. 

lr This formula has already been proposed for generic dimension d, on the grounds of dimensional 
analysis |J. 



6 



One can readily calculate the Id energy 

E d\nZ/L 1 / 8C _\ 1 / dC\ 



L d(3 u /3 2 V 9(3 J (3 2 \ dx 

At zero temperature (/3 — > oo) the system is on its ground state and therefore e 
represents the specific ground state energy whereas the C part is a finite-size effect. 
From the energy, Eq. flT3|), we can compute the thermodynamic entropy 

which is universal, since it contains no contribution from cq. Moreover, the entropy 
vanishes at zero temperature, in accord with the third law of thermodynamics. The 
relation between S and C in Eq. ( fH] ) implies the proportionality between S and c 
at the critical point (CP), namely, S = ttc/(3(3). This is reminiscent of the relation 
between geometric entropy for a CFT and central charge found in |9j . 

The theorem of increase of the relative entropy @ holds on a finite geometry 
and guarantees that S Te i(X, (3, A) increases with A. We calculate the relative entropy 
substituting the value of W = — lnZ/(/?L) = F/L according to Eq. (|l|): 

5rei(A, 0, A) = W(X,0,A) - W(O,0,A) - X dWI - X g ^ A) = 
A > + £ ( C - C < 0) - A f ) = S "'< A ' A) + W - 7< { C - 1 3 • (15) 

where 5 re i(A, A) = lim^oo 5 re i(A, /3, A) is the bulk relative entropy. If this entropy is 
universal we have shown that it takes the form S Te \{X) = B \ 2 / y . Then the presence of 
the scale (3 allows us to define a dimensionless relative entropy, 

C(x) = (3 2 S rel (A, p) = ^ + B x 2 - ( C - - ^ V (16) 

6 \ y ax j 

Furthermore, the monotonicity theorem adopts a dimensionless form, 

dC B 2 f 

x^ = '- d 2 z(:Q(z)::e(0):), (17) 
dx y J 

Since derivatives with respect to x are equivalent to derivatives with respect to (3, C 
embodies RG irreversibility, in the manner of Zamolodchikov's theorem [ rOfl . Although 



C(0) = 0, we can redefine it such that it is proportional to the central charge c at the 
CP by substracting the constant term 7rc/6 from both sides of Eq. flTBD , enhancing 
the similarity with Zamolodchikov's c function. We could say that it also plays the 
role of an off-critical "central charge". From (|16|) it is clear that C(x) has a bulk part 
proportional to x 2 and a finite-size correction, expressed in terms of C(x). As x — » oo, 
C(x) tends to a finite limit and so does the finite-size part of C(x). Hence, in the 
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low-temperature limit x —>■ oo the bulk part dominates, C(x) ~ Bx 2 , so that C(x) 
diverges, unless B = 0. If the relative entropy is not universal, we can nevertheless 
define a dimensionless relative entropy but then as a function of two variables, namely, 
C(x,/3A). Since we must have that (3A ^> 1, monotonicity still holds for moderate 
values of x. 

In parallel with the relative entropy, now it is natural to consider the behaviour of 
the absolute Id quantum entropy S with respect to (3: 

We have again monotonicity, for (3 F is a convex function of f3, as deduced from the ex- 
pression of its second derivative as the average — ((H—(H)) 2 ). Unlike the monotonicity 
of the 2d relative entropy, Eq. (|4]), here H is the total Hamiltonian, that is, including 
the critical part H* — e.g., the kinetic term H* = f(d<p) 2 /2. This monotonicity is in 
principle unrelated to the monotonicity of S rc \ with respect to the coupling constant. 
Thus it allows us to define a different monotonic dimensionless function, 

^/ s S 2C dC 

C{X)= LX^ = -- + ^ < 19 > 

At the critical point S/L = ttc/(3/3), implying that C(x) diverges linearly at x = 0, 
whereas C(0) = 0. On the other hand, as the temperature is lowered (x — > oo) C(x) 
decays to zero. 

We see that there are several quantities that can be related at a RG fixed point but 
have a different physical origin and clearly differ away from it. The quantity which has 
been more prominent in the literature is the finite size correction to the free energy 



C (x) . It was proposed as a monotonic function in Refs. [11], [I| . It has sometimes been 
related to the dimensionless quantity 3 j3 2 (T)/n, which gives the central charge c at 
the fixed point. To clarify this question we prove here that this expectation value is 
instead related to the Id quantum entropy S, showing on the way the general relation 
of expectation values of stress tensor components with thermodynamic quantities. Let 
us consider the expectation values of the complex components of the stress tensor, 
9 := T° and T := Tn — T 2 2 — 2 i Ti 2 , on the cylinder geometry. We have the equalities 

E/L=(T n ), F/L=(T 22 ), 

which come from the definition of the stress tensor and are completely general. One 
deduces that 

S/(Lf3) = (T n -T 22 ) = (T), (20) 

which generalizes the standard relation F/L = (—1/2) (T) [|7], §], actually only valid at 
the fixed point. However, the monotonic function C = (/5/A 1//y ) (T) is related with the 
expectation value (T). In fact, 

m = -l(c-^\ pi) 



(3 2 V 2 dx 
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containing the term xdC/dx, which vanishes at the fixed point. 

The coupling may have been understood in all the above as taking the statistical 
system off criticality. However, nothing in the arguments above requires that for A ^ 
the correlation length be finite. Actually, we can well envisage the situation in which 
a coupling of a system at a multicritical point is such that the coupled system is still 
critical. This situation is described in field theory as a massless flow, which causes 
the system to undergo a crossover ending at another non-trivial fixed point of the RG. 
However, we will only study here massive flows, with a finite correlation length and 
hence a mass parameter m. In free theories, as considered in [[| , m is the mass of the 
particles, bosons or fermions. In interacting theories there is a mass spectrum, which 
can be deduced from the long distance behaviour of the two-point correlation function. 
We will be considering theories soluble with the TBA, which directly renders the mass 
spectrum. One may then select the lowest mass of the TBA spectrum and define the 
dimensionless variable as x = (3 m. In massive theories the function C(x) vanishes 
exponentially as x — > oo, and so do the entropic functions C(x) and C(x). 



3 Finite size thermodynamics for free field models 

3.1 The continuum limit of the lattice Gaussian and Ising 
models 

The 2c? Gaussian model on a square lattice with thermal coupling constant 8 is exactly 
soluble ,Q yielding 

1 f d 2 k 

W($) = - 2 In [1 - 2 6 (cos k x + cos k y )} (22) 

2 J (27r) 

— 7T 

per site |T2|]. It has a CP for B c = 1/4. The continuum limit is performed by redefining 
wave vector as k = ap, a being the lattice spacing, and considering W per unit area. 
Although k belongs to a Brillouin zone, in the continuum limit p runs over the domain, 
—A < PxiPy < A (A ~ vr/a), which becomes the entire plane as A — > oo. In the 
continuum limit we have the field theory of free bosonic particles of mass m such that 

m 2 a 2 = 16 (fi c - B), (23) 

so that y = 2 and the coupling is r = m 2 , omitting an irrelevant proportionality 
constant. 

The relative entropy per unit area of the Gaussian model is best calculated with field 
theory methods, for example, using dimensional regularization PJ . It can be expressed 



2 The thermal coupling constan of a 2d lattice model is of course the inverse 2d temperature. Since 
we shall be using throughout the corresponding Id temperature, (3 — 1/T, we avoid mentioning a 2d 
temperature and use the notation fi for the 2d coupling constant. 



9 



as 

_ r[(4 - d)/2] d/2 

rel ~ (47r)"/2d 

which in d = 2 yields 

5 rc i = (24) 

However, it is more illustrative to start with the expression of the cutoff logarithm of 
the partition function per unit area 

A 

W[r, A] = - In Z[r, A] = \ J -0- 2 In £±I (25) 



which can be integrated exactly and yields 



Naturally, it is UV divergent. For large A it becomes 

W[r, A] = — I -A 2 + r In — + r + 0(A" 2 ) 1 , (27) 
87r [ r J 

exhibiting a quadratic and a logarithmic divergence. 

Recalling the discussion on the general structure of W in the previous section, we see 
that we are in the case of logarithmic corrections to a pure scaling form. Nevertheless, 
it is easily derived that in the present case all the divergences in A cancel in the relative 
entropy, yielding in the infinite cutoff limit 

dW r 

Srei = W{r) - W - r— = — , (28) 

dr 87T 

in accord with the dimensional regularization result. To be precise, in this cutoff 
regularization the quadratic divergence cancels by the substraction of W(0) and the 
logarithmic divergence by the Legendre transform, while in dimensional (or analytic) 
regularization the quadratic divergence does not appear but there is a pole in W, 
equivalent to the logarithmic term in A, that cancels in S Te \. 

Another interesting and exactly soluble example is the 2d Ising model on a square 
lattice, with 

I f d 2 k 

W(S) = -- J -—2 In [cosh 2 (2 6) - sinh(2 fi)(cos k x + cosk y )] - In 2 (29) 



per lattice site [^]. The critical point occurs for the value of B such that the argument 
of the logarithm vanishes when k x ,k y — > 0, namely, when 

/(B) := cosh 2 (2 B) - 2 sinh(2 B) = 0, (30) 



10 



with solution B c = ArcSinh(l)/2 = ln(V2 + l)/2 « 0.440687. The expansion of /(B) 
near fi c yields 

/(fi) = 8 (6 - fi c ) 2 + 8 (fi - fi c ) 3 + 0(6 - fi c ) 4 . 

If we define m by 

m 2 a 2 = 16(6- fi c ) 2 (31) 

and redefine the momentum /c as k = ap, with a the lattice spacing, we obtain near 
the critical point that 

n/a 

W(fl) = ~ / |^ln[(p 2 + m 2 )a 2 /2]. (32) 

— n/a 

In other words, the corresponding field theory is described by a W per unit area given 
by minus that in Eq. (|25[) . It agrees with the well known description of this model in 
terms of a free Majorana fermion theory. However, the relative entropy is not minus 
that of the Gaussian model, since now the coupling constant is proportional to m 
instead of being r = m 2 , since Eq. (|31|) implies that y = 1. One obtains 

S iel (r) = W(r) - W(0) - m^- = (1 + In — ). (33) 

It diverges in the limit of infinite cutoff, which cannot be removed to obtain a universal 
value. Nevertheless, for m <C A, where the field theory makes sense, the relative 
entropy in (|33| ) is monotonic. This is not surprising because it coincides near the CP 
with the exact relative entropy of the square-lattice Ising model, represented in [l|]. 



3.2 Derivation of finite-size quantities 

The expression of W on a lattice of finite size Li x L 2 is obtained by replacing the 
integrals in (^) or ( |2"PD with sums over discrete momenta with step In jh\ and 2n / L 2 . 
When Li,L 2 >awe approach the thermodynamic limit and the sums become integrals 
plus some finite-size corrections. However, the double limit L 1; L 2 — ► oo is complicated 
to study, and it is better to consider finite-size effects only in one direction. Alterna- 
tively, it is sometimes convenient to consider a non-symmetrical lattice with different 
coupling constants in the horizontal and vertical directions. In particular, the quan- 
tum Id Gaussian or Ising models on a chain of sites can be obtained as the extreme 



anisotropic limit of the 2d Gaussian or Ising models [DJ. The CP is still where the 
correlation length diverges but now correlations are calculated only between horizontal 
spins. Now the partition function is Z = Tre _/3H , which can be represented in the 
continuum limit as a functional integral on IR x S 1 with (3 = 1/T the length of the 
compact dimension. It may be good to recall that here (3 has no relation with the 
coupling constant, unlike 8 in the classical 2d models above, and plays instead the role 
of RG parameter. 
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Let us first consider the specific ground-state energy of the Id lattice system. For 
the Gaussian model in the continuum limit it is given by 



A 



/dp X 
— a/p 2 + m 2 = — A V A 2 + m 2 — m 2 log m + m 2 log (A + VA 2 + m 2 ) 
Z7T 47T L 



(34) 

When A — > oo the leading terms are 

e = — (2A 2 + 2m 2 ln— + m 2 + 0(A- 2 )l. (35) 
8ti { m J 

It is quadratically divergent. The logarithmic divergence is universal, that is, indepen- 
dent of the regularization method, and corresponds to the logarithmic divergence of 
VK[r], Eq. (|27l). However, note that the A 2 and m 2 terms are non universal and their 
coefficients change from — 1 and 1 in ( pTf) to 2 and 2 In 2 + 1 , respectively. We shall 
show below that the specific ground-state energy of the Ising model is given by the 
same formula, except for an overall minus sign, in agreement with its free energy in 
Eq. (H). 

In order to compute finite-size effects we first consider the behaviour of the ground 
state energy on a segment of length L at zero temperature, connected with the well- 
known Casimir effect. It provides the finite-size correction C(x) that we need. To see 
this, let us take the specific W, according to Eq. fljjD, 



-InZ C((3m) 

^=eo(A,-) + -^, (36) 

and interchange the roles of L and j3: we have that at low-temperature 

— hiZ , C(Lm) /0 _. 

-e (A,m)+ V r9 ; - (37) 



L/3 uv ' ' L 2 

Since E = —d\nZ/dp, this formula also gives the specific ground-state energy on a 
segment of length L at zero temperature. The term proportional to L is the bulk 
ground-state energy considered above and the finite-size correction is the Casimir en- 
ergy. In other words, the Casimir energy provides the universal function C{mf3). 

3.2.1 Direct calculation of the Casimir energy 

The Gaussian model with periodic boundary condition has a ground state energy 




When L — > oo one recovers the continuum integral of Eq. (0). However, if we are 
interested in the vicinity of the critical theory we may consider the limit L — > oo but 
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with m L small. This limit is known to provide a method to calculate the CFT central 
charge. Then the series can be evaluated by expanding the square root in powers of 
m L and interchanging the sums. We obtain 



m 2-7T 
2 + T 



En 



c(-i; 



m 
~2~ 



1 

+ 2 



2tt 
~L 



z=o 

2 



1/2 



C(i) 



raL 
"2T 

m L 



2/ 



1 



8 V 2vr 



C(2i 



C(3) + 



(39) 



Since £7o is divergent, the result amounts to a zeta-function regularization of it. The 
first term, with ((— 1) = —1/12, yields c — 1. The next term, proportional to L, 
accounts for the bulk term en. Despite the regularization, it is still divergent, since 
((z) has a simple pole at z = 1. This pole is equivalent to a logarithmic divergence 
in regularizations with a UV cutoff, as generally happens when comparing analytic 
with cutoff regularizations. The way to realize it for this case is to restrict the sum 
((1) = YlT(^/ n ) U P t° some large number N. Then 
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(40) 



The connection with the regularization provided by considering the system on a discrete 
chain of spacing a can be made taking iV = L/a, the number of sites. An alternative 
procedure of regularization is first to segregate the divergent bulk part, with the form 
(|34]), from the finite-size corrections by using the Euler-MacLaurin formula (appendix 
A). 

The Ising model on a closed chain is amenable to an analogous treatment. Its ground 
state energy for T > T c is like (EJ3) but with negative sign and with wavenumbers that 



are odd powers of n/L [ITB] 
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The expansion in powers of mL yields 
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The central charge is c = 1/2 and the bulk term is minus that of the Gaussian model. 
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3.2.2 Thermodynamic calculation of finite-size effects 



Now we concern ourselves with the deviation of the energy at non-zero temperature 
from the ground-state energy or, in other words, the finite-size (3 correction to the free 
energy. For the Gaussian model it can be expressed as the free energy of the ideal Bose 
gas constituted by the elementary excitations, 

oo 

Pj = e p+ I gln(l-e-^W), (43) 



where the one-particle energy is e(p) = ^/ p 2 + m 2 . This formula can also be obtained 



by an explicit calculation of the finite-size corrections in the 2c? lattice model |T? . 
When m = it can be used to calculate the central charge Nevertheless, an 
expansion in powers of m? is not advisable: The ensuing integral at the next order is 
IR divergent; that is to say, the expression ( |4"3"D is non analytic at m 2 = 0. Fortunately, 
the integral can be computed by changing the integration variable to e and expanding 
the logarithm in powers of e~ l3e . We obtain 

p 00 i 

P - = e (3 - - Y, - K^nmP), (44) 

71 = 1 

where Ki(x) is a modified Bessel function of the second kind. 

We now define the dimensionless quantity x = m (3 and perform a small-x expansion, 
which yields 

4 = e ° ' - h h -r-T ('■> h + £) + 0(1 ' 4 )} < 45 > 

C(2) m „m 2 /, md l\ ^. 4 . , . 

= ^ - a + 2 + " n 17 + ~< - 2 j + °< m >■ < 46) 

where 0(x 4 ) denotes an analytic remainder of fourth order. The term £(2) = tt 2 /6 
gives the usual m = part and central charge c = 1. The non-analyticity in m 2 of 
the integral for F ( |4*3"D manifests itself in the appearance of the m/2 and logarithmic 
terms. The former also appears as a zero mode contribution in the Casimir energy 
(|38|). The full x-power series is obtained as follows: The series of Bessel functions (|44| ) 
is slowly convergent and one may apply to it a Mellin transform to convert it into a 
rapidly convergent one flBf . Fortunately, its Mellin transform is a power series of x. 
Furthermore, after replacing (3 with L it coincides with Eq. (|39| ) from / = 2 onwards. 
Of course, this should be expected on the grounds of symmetry on a torus under 
interchange of its sides L and f3, that is, modular symmetry. Incidentally, the exact 
free energy on the torus can be computed with some more sophisticated mathematics 



19|, Its two cylinder limits yield Eq. (|39|) or Eq. ([46]). Nevertheless, to show the 



modular invariance of the exact expression on the torus is not easy: It can be done 
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performing its expansion in powers of the dimensionless modular invariant parameter 
m 2 A, with A the area of the torus, but it is very laborious. 

It is interesting to relate the logarithmic term in the expansion of F fl46|) with the 
Casimir energy calculated in the previous subsection. It was remarked there that the 
C(l) divergence can be interpreted as a logarithmic divergence in the cutoff. Adding 
the logarithmic terms in e$ and C{x) one obtains 

x 2 [KB \ x 2 . AT 

where N is the number of lattice sites in the time direction. We see that it is equivalent 
to the divergence ( |4"0"D under the interchange 3 <-> L. 

Now we can calculate the specific entropy 

It has no IR singularity at m — > as opposed to the free energy or the energy. The last 
term is just twice the relative entropy of a box of size L 0, Eq. (p8|), times 0. In more 
generality, for the Gaussian model there is a relation between both types of entropy, 
namely, 

SW(r, P) = W{r, (3) - W(0, d) - r dW ^ P) = S Iel (r) 

+ -[c-C {0 )-r-y-- — + w . (49) 

For other models there is no direct relation between the Id entropy and the 2d relative 
entropy. 

The derivation of thermodynamic quantities for the Ising model is analogous. The 
free energy is that of an ideal Fermi gas 

oo 

/?T = e /3- / ^ln(l + e-^W)=e /5 + -V^K 1 (nm/3) (50) 

^ n=l 
— oo 

where the one-particle spectrum close to the critical point is again e(p) = ^/p 2 + m 2 . 
This integral is computed like the bosonic one. The small-m expansion yields 

C(2) n m 2 / m/3 1\ n , 4 . 

=e ° p -U-'^r^ +i -2) +0{m) - (51) 

As well as for the Gaussian model, it is possible to obtain the x-power series by the 
Mellin transform of the series of Bessel functions (|50|). Similarly, after replacing with 
L it coincides with Eq. (|2|) from I = 2 onwards. 
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In this case the specific entropy is 

For the Ising model the relative entropy is related to the energy, instead of S: 

S, a (m,P) = W(m,fi) - W(0,/3) - m 9W ^ P) = 
*-(") + ^(o-C(0)-mg 

-8^' 1+hl A?»-<L- e »» + l2^- < 53 » 

We see that for free theories we can derive explicit formulas for the free energy— 
and hence for the entropy, — as well as perturbative expansions. Moreover, both the 
Id entropy and the 2d relative entropy give rise to monotonic central charges, as we 
proceed to study, introducing before for convenience the stress tensor. 



3.3 Expectation values of the stress tensor and entropic C 
theorems 

The previous section has shown the calculation of finite-size corrections for various 
quantities of free models with concepts pertaining to the Id quantum theories, namely, 
the lattice Casimir energy or the statistics of quantum gases. The same results can 
be attained with the use of 2d Green function techniques, through the calculation of 
expectation values of the complex components of the stress tensor, taking into account 
their relation with thermodinamic quantities already remarked. We shall rewrite the 
monotonic functions in a suitable way to confirm these relations, hence explaining the 
structure of those functions. We thus start with the expressions of the expectation 
values, 6 := T" and T := Tn — T 2 2 — 2iT i2 , in the cylinder geometry, as derived by 
2d Green function techniques (appendix B): 

m 



(0) = ±— ^o(O) +2j2(±rK (nm(3)j , (54) 

2 / 00 \ 



n=l 



with the same sign convention as before. The modified Bessel functions are divergent 
at zero, namely, K (0) is logarithmic divergent and ^(0) is quadratically divergent. 
These are UV divergences, like those already considered for W, which can be removed 
by normal order. 
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Using the recursion relations satisfied by the Bessel functions we can write the free 
energy (H4T) or fl50D as 



Pj = eo/?T^^(±) n [^ 2 (^m/?)-ir (nm/3)] 

n=l 

= -§ (T-0>=/3(T 22 >, (56) 

showing its relation with the expectation values of the components of the stress tensor, 
an example of the relations obtained at the end of section 1. Notice that it implies a 
definite form for e , to be compared with (|27j) or (|35|). (See appendix B.) 

Similarly, we can calculate 



dW 



9r ^sDw-fl^i 8 '' (57) 



n=l 



771 2 00 

1 = e ° ± ^ E(±) n [^(nm/3) + #0(71 m/3)] 

n=l 

= I(T + e> = (T u >. (58) 

The first equation is just a particular case of the expression of the derivative of W with 
respect to r as the expectation value of the "crossover part" of the action since 9 
is proportional to it. Having the values of W and its derivative available we further 
obtain for the Gaussian model that 

00 

7r r 



5 rel (r,/3) = 5 rel (r) + — - — ^# 2 (nm/3) 



n=l 



5 rc i(r) + ^-i(:T:) (59) 



and for the Ising model that 



00 

7T r v — r 



S ml (r,(3) = 5 rel (r) + — + — ^(-) n [*T 2 (nm/3) + K (nm(3)\ 

n=l 

= ^ei(r) + T ^-i(:T+0:). (60) 

We have substituted S re \{r) for e (r) — e (0) — rde (r)/dr. One obtains the finite 
expectation values of normal-ordered stress tensor components owing to the subtrac- 
tion of W(0, /?).[] Using the connection between the Id entropy and the stress tensor, 
S/(Lf3) = (:T:), pointed out at the end of section 2, one can directly obtain S. 

3 One must be careful when evaluating 9b q^ . Since eo(r) = [-^2(0) — K o (0)] (appendix B), it 

may seem that eo(r) — r de ^ r 7 - = 0. However, -^2(0) and Kq(0) contain an m dependence, because of 
regularization. 
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Thus the dimensionless relative entropies for the Gaussian or Ising models, respec- 
tively, are 



2 2 00 

7T X X 



C(x) = - + > K 2 (nx), 

6 8tt 2tt^ V ; ' 



n=l 
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y J 12 8vr 
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(61) 

f^(-r[if 2 (nx) + K («x)]. (62) 

n=l 

The other monotonic quantity, C = S/(Lm), is common and can be written as 

00 

C(x) = -Y(±) n+1 K 2 (nx), (63) 

IT < ^ 



n=l 



which is, of course, C = ((3/m) (:T:), according to the expression of (T) (p>5|). Series 
expansions of C, C and C are derived from Eq. (|6|) or (|5T|). Both C and C are plotted 
in Fig. 1. The Ising model C is for the value (A f3) 2 = 10000. It is useful to recall that 
in general (A/?) 2 ~ N, the number of 2d lattice sites in a box of side (3. 
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Ising model 
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Fig. 1: C and C for the Gaussian and Ising models. In the latter model C is for A/3 = 100. 



4 Interacting models 



For interacting models, the free energy is in principle not available in closed form. 
Nevertheless, one can perform its perturbative expansion. In two dimensions one can 
take advantage of the information provided by the methods of conformal field theory 
(CFT), namely, the non-perturbative dimensions and correlations of fields at the critical 
point; then one speaks of deformed CFT's. For example, one can perform a perturbative 
expansion around the critical point. This approach, called conformal perturbation 
theory, is well suited for the calculation of an entropy relative to the critical point; 
the perturbation parameter is A (3 y . The expansions of Eq. ( ]46|) and ( fTl]) are instances 
of it and can be obtained from the respective c = 1 or c = 1/2 CFT [15]. The 



logarithmic term in them is due to UV divergencies. The analysis of UV divergences 
is done by examining the behaviour of integrals of correlators for coincident points. If 
1 < y < 2 there are no UV divergences in the perturbative expansion and, furthermore, 
this expansion is arguably convergent ||, ||. in contrast, when y < 1 a finite number of 
terms will diverge. The condition 1 < y < 2 agrees with the non-perturbative regularity 
condition for the relative entropy found before. Conformal perturbation theory is 
considerably powerful but, at any rate, the perturbative expansion only converges for 
a limited range of A /3 y , while we are interested in the behaviour of thermodynamical 
quantities over the entire range of the coupling constant. 

Some 2d models are partially soluble with the thermodynamic Bethe ansatz (TBA) 
20fl . In particular, models for which the interaction is of purely statistical nature lend 



themselves to a derivation of closed expressions for the free energy and entropy similar 
to the ones for free field models, albeit more complicated. Hence, the entropic C the- 
orems can be explicitly verified for them. The Bethe ansatz assumes a factorized form 
for the wave functions and hence an expression for the energy as a sum of contributions 
of independent quasi-particle levels, though these quasi-particles have non-trivial ex- 
change properties. To determine the structure of these levels is a complicated business 
but it dramatically simplifies in the thermodynamic limit, constituting the basis of 
the TBA. This method is non perturbative in nature and provides thermodynamical 
quantities over the entire range of the coupling constant. In principle, the TBA yields, 
further to the universal finite-size correction function C(x), a contribution proportional 
to x 2 , which is interpreted as a UV-finite bulk term and therefore has been called uni- 
versal bulk term |21], Therefore, in this section we redefine C(x) to include this 
universal bulk term. 

In spite of the virtues of the TBA approach, the TBA equations themselves are by 
no means easy to solve and it is customary to resort to numerical calculation to obtain 
the coefficients of the series expansion in (5 m. In this sense the TBA approach is not 
superior to perturbation theory in its region of convergence, with which one obtains 
analytic expressions for these coefficients. It is only in the case of purely statistical 
interaction where the TBA approach is definitely superior, for one can then solve 
the TBA equations algebraically. Thus we treat this case first. It applies to models 
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of Calogero-Sutherland type, which represent the dynamics of spinons or other non- 
interacting particles with fractional statistics. We can calculate the entropy and check 
that is monotonic with respect to (3. However, it is beyond our means to calculate 
expectation values of the stress tensor, since the expression of the stress tensor is not 
available. In second place, we shall treat the general interacting integrable case by the 
numerical solution of the TBA equations and compare with the results of conformal 
perturbation theory. 

Last, there is an alternative to conformal perturbation theory or the TBA that can 



be applied to any model, namely, numerical finite-size scaling on a chain ||22|| . Like 
the TBA, by its own nature it is not limited to a restricted range of (3 m. However, 
the numerical calculations required to obtain similar accuracy to that of the TBA 
are prohibitive in practise. Thus this brute-force method is not actually effective to 
compute off-critical quantities and no use will be made of it here. 



4.1 Models with purely statistical interaction 

The role of fractional statistics in Condensed Matter Physics, as a generalization of 
the regular bosonic or fermionic symmetry properties under particle exchange, has 
been long recognized [^]. Its modern version has given rise to the concept of anyons. 
As it happens, this type of statistics leads to highly non-trivial correlations between 
particles which are difficult to disentangle and indeed constitute what has been called 
statistical interaction. The form of this interaction can be best realized by transforming 
the particles to standard fermions or bosons with a peculiar interaction. Models with 



purely statistical interaction are usually referred to as generalized ideal gases |24|, |25 
It is customary to consider the free particles as fermions and parametrize the statistics 
by a number g, such that the maximum number of particles that can fit in a single 
fermion momemtum level is 1/g. For no statistical interaction g — 1. If g — 1/n 
the single-fermion levels can accomodate n particles and in the limit n — > oo the 
statistics becomes bosonic. Some models with apparently complex interactions can be 
transformed into generalized ideal gases, as occurs for the Calogero-Sutherland models 
26] or their lattice version J27] . 



Therefore, models with purely statistical interaction are interesting systems and, in 
addition, sufficiently complex to be a suitable benchmark for our irreversibility theorem. 
We start with the spin SU(2) level one Wess-Zumino-Witten Lagrangian in the bosonic 
representation, 

H = j dx[{d t cb) 2 + {d x cb)\ (64) 
The spinon field is defined in terms of the bosonic field as tp — ^P^^f wnere 



the sign stands for spinon polarization f27| , |28| , p9| , |30| . It is a free theory, which 
can also be expressed as a free fermionic theory, but we are going to do some non- 
trivial manipulations on it. First, we "simplify" the model by keeping just one spinon 
polarization, "+" say. The physical way to achieve this is to introduce a very strong 
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magnetic field. Now we have a semionic CFT with central charge c = 3/5 |3T|, |32 



which is certainly an interacting theory. However, its partition function is known, as the 
total partition function of the SU(2) level one Wess-Zumino-Witten mode restricted 



to vanishing fugacity of "— " spinons, z_ = 0, [29]. The thermodynamic quantities can 
be obtained with the help of the TBA.Q 

The second change consists of the addition of some tunable coupling, which perturbs 
the model away from criticality and allows one to probe the behaviour of the entropy. 
If we impose that the interaction remains purely statistical, the only possibility is to 
give mass to the semions: We can replace the dispersion relation e(p) = \p\ with e(p) = 
\/p 2 + m 2 . We assume that this perturbation fulfills the conditions for the application 
of the monotonicity theorem ([|) and we shall see that it is easily implemented within 
the TBA approach and yields expressions which can be treated by algebraic methods. 
As a side remark, note that the massive relativistic dispersion relation differs from the 
non-relativistic one assumed in the original Calogero- Sutherland model. Therefore, if 
we want to keep to the physics represented by this model, we must interpret the mass 
as a parameter unrelated to the real semion mass. For comparison, it can be shown 
that the 2d Dirac Lagrangian with a mass term appears as a low energy effective 
Lagrangian for (non-relativistic) conducting electrons in one-dimensional metals, but 
the Dirac mass is actually related to the electric potential |3(], chapter 13]. Regardless 
of the precise physical interpretation, we will consider the theory of relativistic massive 
semions as our first interacting field theory to investigate the properties of the entropy. 



4.1.1 Application of the TBA to the semion gas 



The full power of the TBA shows in the calculation of finite size corrections to ther- 
modynamic quantities. One obtains for the critical semion gas |32], [33| 
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e(k) = \k\. We give the semions a mass, replacing the dispersion 



relation e(k) = \k\ with e(k) = y/k 2 + m 2 . Then 
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4 Derivation of thermodynamic quantities from CFT usually demands a thermodynamic approach 
in the sense of 32 , be the TBA or Schoutens' recursion method |n) . 
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One can easily obtain its perturbative expansion in powers of m (3 by expanding first 
the modified Bessel function, 
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where f/>(x) is the digamma function. However, we will content ourselves with extract- 
ing the bulk part, 
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ln(2n + l) = -0.0536114. 



Corresponding to g = 1/2 statistics, in a segment of length L the semion momenta 
of completely filled single-particle levels are k n = (2tv/L) |n/2 + l/4| with n G Z. Hence, 
the approximate ground state energy is like ([II]) but with wavenumbers that are odd 
powers of 7r/(2 L) |31], |32|, 
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Now, an expansion in powers of mL would yield a wrong Casimir energy, owing to the 
approximated nature of the ground state energy. We can however see that the bulk 
term 
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is twice that of fermions, as corresponds to the double average occupation number 
of semions. We can compare the approximate bulk non-analytic term in (|7T|) with 
the exact result of the TBA (|69|) . According to the combined form of IR and UV 
logarithmic terms ( f4"T|) that we expect from (|7T|), 



C{x) |buik = In x + 0.390536 x 2 . 
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We see that the first term of the sum Si reproduces the coefficient in this approximation, 
1/2, but the total coefficient, l/VE ~ 0.447214, is slightly smaller. 

The total non-analytic part is now an infinite series, obtained from 
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Interestingly, the series coefficients seem to be rational numbers. 

The entropic C functions are easily derived from the expression of C 
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They are plotted in Fig. 2. We use again the value A/3 = 100. 




Fig. 2: C (A (3 = 100) and C for semions. 



4.2 Deformed two-dimensional conformal field theories 

A general class of theories amenable to derivation of the finite-size quantities of interest 
is that of deformed 2d CFT. Since CFT provides the exact dimensions of relevant fields, 
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the results of conformal perturbation theory are more accurate in 2d than those of 
ordinary perturbation theory, which on the other hand is plagued with IR problems. 
In addition, many models admit integrable deformations, in the sense that the existence 
of an infinite number of conservation laws forces the S-matrix to factorize. Then the 
TBA provides a way to derive thermodynamic quantities. Many deformed CFT are 
known to be integrable and similar methods are applicable to all f34" l, although their 



complexity can be considerable for the most sophisticated models. Therefore, we shall 
choose one of the simplest cases. It should be intuitively clear how to generalize the 
computation of the entropic quantities to other integrable models. 

The natural (and oldest) generalization of the Ising model consists of taking a site 
variable which can take three values instead of two, constituting the three-state Potts 
model. It is critical for B c = ln(\/3 + l)/3 and its thermal critical exponent is v = 5/6, 
implying that y = 6/5. This model has been long known to be integrable and it has 
been long (but not as long) known to be describable in terms of particles with fractional 
statistics, which are generalizations of the Ising fermions and are called parafermions. 
This quasi-fermionic representation is in terms of two conjugate parafermions carrying 
2 3 charges +1 or —1 and spin 2/3, which are massless at the CP point but acquire a 
mass for T > T c . Their interaction is purely statistical at the CP but it is more com- 
plicated off criticality. However, it is still integrable and its thermodynamic properties 
can be found with the TBA. It yields 



C(/3m) = -(3m 



oo 

/rlf) 
— cosh(fl) ln(l + e- e(e) ) , (76) 
27T 



which is apparently similar to the formula for free fermions but now e{6) are unknown 
functions to be determined with the TBA equations. The concrete ^-matrix elements 
of this model lead to the TBA equation 

oo 

<(w°*'+¥ I ^^VV t 1 ^)- (77) 

— oo 

It can be solved numerically by an iterative algorithm, yielding a set of numbers which 
can be displayed in a table pifl . Hence, according to the general formula (|i~6|), 

„, . Ti c . 5x dC(x) . 

C(,) = T+ 0(x)- T - 5 U, (78) 

where the bulk relative entropy S re \(m) does not appear since it is now implicitly 
included. It derives from the universal bulk free energy, which can be calculated exactly, 
yielding -y/3x 2 /6 |2T|; hence, S Te i(m) = V3m 2 /9. We can calculate the monotonic 
functions C C numerically as well. They are plotted in Fig. 3. 
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3-Potts model 
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Fig. 3: C and C for the 3-state Potts model. 



The previous solution for C(x) can be expressed as a power series in x 12 / 5 plus a 



bulk term 21 
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Its radius of convergence can be estimated to Sx 12/5 = 14.3 ± 0.4, that is, Sx ~ 3.0.0 
This expansion ( |79"|) can also be obtained by conformal perturbation theory pT| . The 
series coefficients are then expressed in terms of integrals of correlators of the perturbing 
conformal field, that is, the thermal field of dimension <i$ = 4/5 in the three-state Potts 
model. The computation of these integrals is very laborious, except for the first ones, 
which in some cases admit explicit expressions as series A related computational 
method is the truncated conformal- space approach, in which one only takes a finite 
dimensional subspace of the Hilbert space of possible states, namely, a number of low- 
lying states, to numerically diagonalize the Hamiltonian. In the limit of increasing the 
number of states, this approach is equivalent to a numerical evaluation of the integrals 
giving the coefficients of the perturbation series 0. We do not think worthwile to dwell 
in detailed computational methods since the simple numerical integration of the TBA 
equation Q77D suffices to show the monotonicity of C and C. 



5 A plot of this series sharply shows that it is very close to the TBA result for x < 3 and quickly 
departs from it for larger x. 1 1 
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5 Discussion 



For a 2d field theory, one can introduce two monotonic dimensionless functions, namely, 
C and C, derived from the 2d relative entropy and the Id quantum entropy, respectively. 
It has been shown that C is universal when the couplings are strongly relevant, that 
is, with dimension 1 < y < 2. They include the thermal perturbations of the unitary 
minimal models of conformal field theory (CFT), except for the Ising model, which we 
have also studied, notwithstanding. In contrast, C is always universal, since it only 
depends on the universal finite-size correction to the free energy. Given that general 
theorems may not be particularly useful if the quantities that they involve cannot be 
computed in practise, considerable time and effort has been devoted to compute C and 
C for a variety of models. In consequence, we have been able to show for them that 
those functions are monotonic. 

The dimensionless entropies C and C play a similar role to Zamolodchikov's c func- 
tion, constraining the structure of the RG flow, but they have a clear physical origin, 
unlike Zamolodchikov's c function. The existence of a monotonic function is usually 
argued on the grounds of the irreversible nature of the RG flow, which in the coarse- 
grained formulation implies a loss of information on microscopic degrees of freedom 



T0| , |35| . This idea inspired the adaptation of Boltzmann's H theorem to the RG flow 
in our previous work It has been shown here that this philosophy gives rise to the 
entropic functions C and C, which are computable for a wide range of models. The non- 
perturbative computation of Zamolodchikov's c function is much harder and in fact it 
does not seem to have been carried out for any fully interacting model. For all these 
reasons, the entropic C functions proposed here arguably provide a new perspective in 
the long-standing problem of the irreversibility of the RG. Nevertheless, in compari- 
son with Zamolodchikov's c theorem, it must be remarked that universality of C, the 
entropic function more similar to Zamolodchikov's c function, has been proved only 
for deformations of the critical theory by fields with dimension < d$ < 1 (strongly 
relevant), while Zamolodchikov's theorem covers the entire range of d$. 

Our entropic monotonicity theorem for the dimensionless relative entropy is 

= — Id 2 z(:Q{z)::Q{0):). (80) 
dx y J 

Even though it resembles Zamolodchikov's c theorem it is not quite the same: The 
correlator of G's in the second term appears integrated. Furthermore, a detailed cal- 
culation of Zamolodchikov's function c(m) for the free boson or fermion shows that 
they differ from the respective values of C(m) = C(x)\/3=i. The essential discrepancy 
actually has a geometrical origin: A crucial step in the proof of Zamolodchikov's the- 



orem relies on the assumption of rotation symmetry pl| , which does not exist on the 



cylinder. Therefore, the theorem does not hold on it. However, the absence of rotation 
symmetry is traded for the appearance of a new parameter, the width (3, which replaces 
the distance to the origin in Zamolodchikov's theorem and is used in the derivation of 
the entropic monotonicity theorems. 
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The reason for the introduction of a finite geometry is to have an IR scale to define 
a dimensionless relative entropy. We have used the cylinder because of its thermo- 
dynamic interpretation. Of course, other finite geometries are possible. For example, 
one can use a sphere. Its radius is then the IR scale. The advantage is that rotation 
symmetry is preserved on the sphere and Zamolodchikov's theorem holds. With this 
new geometry the monotonicity theorem for the dimensionless relative entropy would 
still involve an integral of the correlator of 6's but a relation of C with Zamolodchikov's 
c function seems more feasible. At least in conformal perturbation theory one should 
be able to perform that integral in terms of the IR scale and a direct comparison with 
Zamolodchikov's theorem could be possible. 

The existence of several monotonic functions prompts the question of which one is 
preferable. It is intuitively clear that a unique definition of monotonic function is not 
possible: The RG itself is not unique and one can choose a variety of RG parameters. 
Correspondingly, if C(x) is monotonic a monotonic change of the independent variable 
x will transform it into a different monotonic function. We might then consider what 
happens at the boundary, x = or x — > oo. The point x = is the RG fixed point and 
it is sensible to define a function related to it, as are the dimensionless relative entropy 
C or Zamolodchikov's c function. The dimensionless absolute entropy C(x) is defined 
irrespective of the fixed point and actually diverges there. We could as well demand 
good behaviour in the limit x — > oo. This condition is fulfilled by C but may not be 
fulfilled by C, owing to the bulk term. It is quite possible that a minor modification 
of the definition of C may remove the bulk term and make it well behaved in the limit 
x — > oo as well as at x = 0. 
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A Calculation of finite-size corrections with the Euler- 
MacLaurin formula 



For free models the energy series can be evaluated with the Euler-MacLaurin summa- 
tion formula,^ 

oo 

t «W = / *■ «W +5 ««» " T2 f '<°» + 4 £ "' (0) " 30^40 £< " ,(0) + " ' ' (81) 

n=0 J Q 

where e(n) are single-particle energies. For the Gaussian model the first term can be 
proved to be proportional to L with the change of variable p = 2tt n/L and leads to the 
integral in Eq. fl34|) . One can see that all the odd derivatives of e(n) vanish at n = 
because it is an even function. It is natural, because the finite size corrections are 
exponentially negligible when L — > oo and hence non analytic. Since every derivative 
pulls out a power of 2n/L, the subsequent series of powers oil/ L must have vanishing 
coefficients. 

In the scaling zone mL < 1 the finite-size corrections are not negligible, despite 
the previous argument. To evaluate these corrections we can nevertheless make use of 
the Euler-MacLaurin expansion but using a non-zero value for the point at which the 
derivatives are computed, in the following form, 

oo 00 

yVn)= [ dne(n) + ~e(l)-— e'(l) + — e"'(l) — e {v) (1) + ■ ■ ■ . (82) 

v ) j v j 2 v ; 12 v ; 72Q u 30240 v ; v ) 

n=l 

Now the series can be cast as an expansion in powers of mL. The reason why this 
trick works can be understood in several ways. One is that the derivatives e^ 2k+1 \n) 
as functions of m are ill behaved for small n. They converge to the null function for 
n = but non uniformly. It is actually safer to choose the argument n of the derivatives 
larger than in Eq. fl82|) , n = 3 or 4 say. Then the Euler-MacLaurin expansion converges 
very fast and the terms displayed above suffice to match the coefficients in Eq. (|39|) 
with about ten decimal places. 

Alternatively, one may focus on the fact that for m = the function e(n) oc |n| is 
singular at n = 0; its derivatives eventually diverging there. This is naturally an IR 
divergence, which does not exist for m^O. However, one must be careful to evaluate 
e(n) at n ^ before taking m = 0, or in other words, one must introduce an IR cutoff 
and evaluate the Euler-MacLaurin expansion at that point. Of course, the result shall 
be inpependent of the precise value of the cutoff, although its convergence properties 

6 This is a common method to convert sums to integrals. However, since the function summed e(n) 
diverges when n — > oo a preliminary regularization is required. A convenient form is to sum up to 
some arbitrary number N 3> 1, which for a chain can be the number of sites. This UV regularization 
renders meaningful the formal manipulations that follow. However, we do not need to be definite on 
the UV regularization for our focus is on universal quantities. 
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are greatly affected by it. Within the realm of classical mathematics, it is interesting 
to recall that Legendre met a similar problem when he attempted to evaluate elliptic 
integrals numerically with the Euler-MacLaurin expansion. Since the integrand is an 
even function of the integration variable at the limits and tt/2, the odd derivatives 
vanish and the Euler-MacLaurin formula implies that the elliptic integral is equal to any 
of its rectangular approximations. The paradox was solved by Poisson, who showed that 
in this case the remainder term does not tend to zero as the number of terms increases 
and hence the series does not converge. If we further consider that the bulk free energy 
of the Gaussian or Ising models on a finite chain can be expressed as elliptic integrals, 
we may appreciate that Legendre actually encountered an IR divergence without being 
aware of the need of regularization. 

For the Ising model fl4"T| ) the odd derivatives of e(n) do not vanish at n = but it 
also is necessary to choose n = 3 or 4 for fast convergence. 



B Calculation of (Tab) on the cylinder for free mod- 
els 



For free field theories the expectation values of the components of the stress tensor can 
be expressed in terms of the Green function. Thus for a bosonic field 

(6) = m 2 (<p 2 ) = m 2 lim G (z, z). (83) 

2— >0 

We use complex notation, z = X\ + ix 2 - The Green function on a cylinder, Gp(z, z) is 
non trivial. Its Fourier transform includes a sum over discrete momenta in the compact 
direction, 

oo 

1 . . r rlh pi (^n xi+k x 2 ) 

a / — ' / 27T uit + k z + m l 

n=— oo " 
— oo 

where the allowed frequencies for bosons are u n = ^ n. It can be transformed into 
a more manageable form by the use of the proper-time representation |36|]. In this 
representation the integral over k is elementary and one is left with the sum over n and 
the integral over proper time. After performing a convenient Poisson resummation one 
obtains 

oo 00 

Ef ds 2 |z-«/3| 2 

hrs e ~ m (85) 

n=— oo q 

Since the Green function on the plane is just 

G^z^z) = —K (m\z\), 

Z7T 



which is the n = term in the sum (p5|), this sum can be interpreted as the solution 
of the field equation for a point source by the method of images. It can be expressed 
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in terms of the Jacobi theta function 

00 2 



n=— 00 



as 

00 



Airs \ Airs Airs 







Then the Poisson resummation realizes the duality property of 63. 
The formal z — > limit of Gp is easily taken, 



00 „ 

G p (0) = ^mG p (z,z)= J2 J 



ds m 2„ (n^) 2 

e 4 * . 



47TS 



n=— oo 



It contains the logarithmic divergence 



00 



G„(0) = ^„(0)=/A e -~\ 







Taking into account the integral representation of modified Bessel functions 



00 







one obtains 

j 00 

^(0) = — E Kv(nmP). 



2tt 

n=— 00 



The computation of (T) requires a little more work, for 

(T) = -4((^) 2 )=4 1im^). 

z^O 



From 

d 2 z G,(z,z)= J2 

Hence, 



n=— oo 



/ \ 2 

ds fz + n[3\ 2 \z-nf « 2 



Airs \ As 



OO 

ajG(0) = limSjG W = V f^^l e -^-^ = -'f ± K A r, 

Z Z^O 2 ^ / 47TS 16S 2 87T ^ 

n=— 00 p, n=— 00 
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It contains a quadratic divergence in K 2 (0). 

For Majorana fermions one could start from the known expressions of their stress 
tensor but it is simpler to consider them as bosons with antiperiodic boundary condi- 
tions and use again the Fourier transform (|4j) with allowed frequencies u> n = ^ (n + |) . 
Before Poisson resummation Gp can be expressed in terms of 

oo 
n=— oo 

The Poisson resummation transforms 9 2 into its dual 84, 

00 2 
4 (i/, r)= (-)"«"?*, 



which is like the bosonic #3 but with an additional alternating sign. 
We can write the final result in a condensed notation, 



9 00 

m 



<e> = ±— J2 (±) n K (nmP), (93) 



n=— 00 
2 00 



( T ) = ± V~ E (±) n ^(nm/?), (94) 

n=— 00 

where the upper signs stand for bosons and the lower signs for fermions. 

However, to have well defined expressions we must further introduce a regularization 
that removes the divergences in K (0) and K 2 (0). It is customary to begin defining 
normal-ordered composite fields, namely, 

:0: = m 2 :^ 2 :, (95) 
:T: = -4:(«9^) 2 :, (96) 

in the sense of a point splitting regularization and a substraction of the divergent part, 
computed with the Wick prescription. It amounts to the substraction of Goo(0) = 
±1/(2tt) K o (0) or 4<9 2 Goo(0) = ±m 2 /(2n) K 2 (0). Point splitting on a lattice yields 

ir (0) = -(ln^ + 7 )+O(a 2 ), (97) 

K 2 {U) = ^-\ + 0{a 2 ). (98) 

Considering a ~ 1/A we have, for example, that according to Eq. fl56| ) 

e = ±^ [JTo(O) - K 2 (0)} ~ ±-L ( -2 A 2 + m 2 In — + m 2 ( 7 - ^ + 0(A" 2 ) 

(99) 



31 



References 

[1] J. Gaite and D. O'Connor, Field Theory entropy, the H -theorem and the renor- 
malization group, Phys. Rev. D54 (1996) 5163; J. Gaite, seminar at the RG96 
Conference, Dubna, [hep-th/9610040| 

[2] J. Gaite, Relative entropy in 2d QFT, finite-size corrections and irreversibility of 
the RG, Phys. Rev. Lett. 81, 3587 (1998) 

[3] A.H. Castro Neto and E. Fradkin, The Thermodynamics of Quantum Systems and 
Generalizations of Zamolodchikov's C-theorem, Nucl. Phys. B400 (1993) 525 

[4] R. Guida and N. Magnoli, All order I.R. finite expansion for short distance behav- 
ior of massless theories perturbed by a relevant operator, Nucl. Phys. B471 (1996) 
361-388 

[5] T.R. Klassen and E. Melzer, The thermodynamics of purely elastic scatter- 
ing theories and conformal perturbation theory, Nucl. Phys. B350 (1991) 635; 
T.R. Klassen and E. Melzer, Spectral flow between conformal field theories in 1 + 1 
dimensions, Nucl. Phys. B370 (1992) 511 

[6] F. Constantinescu and R. Flume, The convergence of strongly relevant perturba- 
tions of d = 2 conformal field theories, Phys. Lett. B326 (1994) 101-104 

[7] H.W.J. Blote, J.L. Cardy and M.P. Nightingale, Conformal invariance, the central 
charge and universal finite-size amplitudes at criticality, Phys. Rev. Lett. 56 (1986) 
742-745 

[8] I. Affleck, Universal term in the free energy at a critical point and the conformal 
anomaly, Phys. Rev. Lett. 56 (1986) 746-748 

[9] C. Holzhey, F. Larsen and F. Wilzek, Geometric and renormalized entropy in 
conformal field theory, Nucl. Phys. B424 (1994) 443-467 

[10] A.B. Zamolodchikov, Irreversibility of the flux of the renormalisation group in a 
2d field theory, JETP Lett. 43 (1986) 730-732 

[11] D. Boyanovsky and R. Holman, The c-theorem revisited, Phys. Rev. D40 (1989) 
1964 

[12] G. Parisi, Statistical Field Theory, Addison- Wesley Pub. Co., 1988 

[13] L. Onsager, Crystal statistics I: A two-dimensional model with an order disorder 
transition, Phys. Rev. 65 (1944) 117 

[14] M. Suzuki, Equivalence of the two-dimensional Ising model to the ground state of 
the linear XY model, Phys. Lett. A34 (1971) 94 



32 



[15] E. Fradkin and L. Susskind, Order and disorder in gauge systems and magnets, 
Phys. Rev. D 17 (1978) 2637 

[16] T.D. Schultz, D.C. Mattis and E. Lieb, Two-dimensional Ising model as a soluble 
problem of many fermions, Rev. Mod. Phys. 36 (1964) 856 

[17] C. Nash and D. O'Connor, Modular invariance of finite size corrections and a 
vortex critical phase, Phys. Rev. Lett. 76 (1996) 1196-1199 

[18] G.G. MacFarlane, The application of Mellin transforms to the summation of slowly 
convergent series, Philos. Mag. 40 (1949) 188 

[19] H. Saleur and C. Itzykson, J. Stat. Phys. 48 (1987) 449; C. Itzykson and J. Drouffe, 
Statistical field theory, vol. II, Cambridge University Press, (1989). 

[20] C.N. Yang and CP Yang, Phys. Rev. 147 (1966) 303; J. Math. Phys. 10 (1969) 
1115 

[21] ALB. Zamolodchikov, Thermodynamic Bethe ansatz in relativistic models. Scaling 
three state Potts model and Lee-Yang models, Nucl. Phys. B342 (1990) 695 

[22] P. Christe and M. Henkel, Introduction to conformal invariance and its applica- 
tions to critical phenomena, Springer Lecture Notes in Physics, Vol. ml6 (1993) 

[23] E. Fradkin and L.P. Kadanoff, Disorder variables and parafermions in two- 
dimensional statistical mechanics, Nucl. Phys. B170 (1980) 1-15 

[24] Y.-S. Wu, Statistical Distribution for Generalized Ideal Gas of Fractional- Statistics 
Particles, Phys. Rev. Lett. 73 (1994) 922 

[25] S.B. Isakov and S. Viefers, Model of statistically coupled chiral fields on the circle, 
Int. J. Mod. Phys. A12 (1997) 1895-1914 

[26] F. Calogero, Solution of a three body problem in one- dimension, J. Math. Phys. 10 
(1969) 2191-2196; Ground state of one- dimensional N body system, J. Math. 
Phys. 10 (1969) 2197-2200; 

B. Sutherland, Quantum many-body problem in one- dimension: ground state, J. 
Math. Phys. 12 (1971) 251 

[27] F.D.M. Haldane, Exact Jastrow-Gutzwiller resonating-valence-bond ground state 
of the spin-(l/2) antiferromagnetic Heisenberg chain with 1/r 2 exchange, Phys. 
Rev. Lett. 60 (1988) 635; 

B.S. Shastry, Exact solution of an S—l/2 Heisenberg antiferromagnetic chain with 
long-ranged interactions, Phys. Rev. Lett. 60 (1988) 639 

[28] P. Bouwknegt, A. Ludwig and K. Schoutens, Spinon Bases, Yangian Symmetry 
and Fermionic Representations of Virasoro Characters in Conformal Field Theory, 
Phys. Lett. B338 (1994) 448-456 



33 



[29] F.D.M. Haldane, Z.N.C. Ha, J.C. Talstra, D. Bernard and V. Pasquier, Yangian 
symmetry of integrable quantum chains with long-range interactions and a new 
description of states in conformal field theory, Phys. Rev. Lett. 69 (1992) 2021; 
D. Bernard, V. Pasquier and D. Serban, Spinons in Conformal Field Theory, Nucl. 
Phys. B428 (1994) 612-628 

[30] A.M. Tsvelik, Quantum Field Theory in Condensed matter Physics, Cambridge 
University Press, (1995). 

[31] K. Schoutens, Exclusion Statistics in Conformal Field Theory Spectra, Phys. Rev. 
Lett. 79 (1997) 2608-2611 

[32] J. Gaite, On thermodynamic approaches to conformal field theory, Nucl. Phys. 
B525 (1998) 627-640 

[33] A.G. Bytsko and A. Fring, Thermodynamic Bethe ansatz with Haldane statistics, 
Nucl. Phys. B532 (1998) 588-608 

[34] G. Mussardo, Off-critical statistical models: Factorized scattering theories and 
bootstrap program, Phys. Rep. 218 (1992) 215-382 

[35] J.L. Cardy, Conformal invariance and Statistical Mechanics, in Fields, strings and 
critical phenomena, Les Houches 1988, Ed. E. Brezin and J. Zinn- Justin, Noth- 
Holland, Amsterdam, 1990 

[36] H.W. Braden, Expansions for field theories on S 1 x E, Phys. Rev. D25 (1982) 
1028 



34 



